Preprint typeset using WT^^ style emulateapj v. 8/13/10 



A STATISTICAL STUDY OF THE LUMINOSITY GAP IN GALAXY GROUPS 

Saeed Tavasoli 

Department of Physics, Ferdowsi University of Mashhad, Mashhad, Iran and 
School of Astronomy and Astrophysics, Institute for Research in Fundamental Sciences (IPM), P. O. Box 19395-5531, Tehran, Iran 



Habib G. Khosroshahi 

School of Astronomy and Astrophysics, Institute for Research in Fundamental Sciences (IPM), P. O. Box 19395-5531, Tehran, Iran 

Ali Koohpaee 

School of Astronomy and Astrophysics, Institute for Research in Fundamental Sciences (IPM), P. O. Box 19395-5531, Tehran, Iran and 
Faculty of Science, Khajeh Nasir Toosi University of Technology, Tehran, Iran 

Hadi Rahmani 

Inter-University Center for Astronomy and Astrophysics, Post-Bag 4, Ganeshkhind, Pune 411007, India 

AND 

Jamshid Ghanbari 

Department of Physics, Ferdowsi University of Mashhad, Mashhad, Iran 

ABSTRACT 

The luminosity gap between the two brightest members of galaxy groups and clusters is thought 
to offer a strong test for the models of galaxy formation and evolution. This study focuses on the 
statistics of the luminosity gap in galaxy groups, in particular fossil groups, e.g. large luminosity 
gap, in an analogy with the same in a cosmological simulation. We use spectroscopic legacy data of 
seventh data release (DR7) of SDSS, to extract a volume limited sample of galaxy groups utilizing 
modified friends-of-friends (mFoF) algorithm. Attention is paid to galaxy groups with the brightest 
group galaxy (BGG) more luminous than Mr— -22. An initial sample of 620 groups in which 109 
optical fossil groups, where the luminosity gap exceeds 2 magnitude, were identified. We compare the 
statistics of the luminosity gap in galaxy groups at low mass range from the SDSS with the same in 
the Millennium simulations where galaxies are modeled semi-analytically. We show that the BGGs 
residing in galaxy groups with large luminosity gap, i.e. fossil groups, are on average brighter and live 
in lower mass halos with respect to their counter parts in non-fossil systems. Although low mass galaxy 
groups are thought to have recently formed, we show that in galaxy groups with 15 galaxies brighter 
than Mr > —19.5, evolutionary process are most likely to be responsible for the large luminosity gap. 
We also examine a new probe of finding fossil group, Amu > 2.5, and find that the fossil group 
selected according to new probe are more abundant than those selected using the conventional probe, 
Ami2 > 2, in low halo mass regime, < IQI^Mq . In addition we extend the recently introduced 
observational probe based on the luminosity gap, the butterfly diagram, to galaxy groups and study 
the probe as a function of halo mass. This probe can, in conjunction with the luminosity function, 
help to fine tune the semi-analytic models of galaxies employed in the cosmological simulations. 
Subject headings: galaxy groups: fossil, nonfossil, richness, galaxy: luminosity function, luminosity 
gap 



1. INTRODUCTION 

Galaxy groups are key systems in advancing our un- 
derstanding of structure formation and evolution in the 
universe as, in a hierarchical framework, they span the 
regime between individual galaxies and massive clus- 
ters. Advances in the cold dark matter cosmological 
simulations and improvements on the particle mass res- 
olution has provided the opportunity to study galaxy 
groups or low mass halos. The Millennium simulation 
(jSpringel et al.l 120051 ) is a recent example. The semi- 
analytic models of galaxy formation are also employed 
to understand the formation and evolution of galaxies in 
cosmological context. 

|tavasoli@ipm.ir] 



The key factor in the success of this attempt is the 
availability of observational constraints to help fine tune 
the models. While galaxy luminosity function has been 
widely used as a strong constraint, the number of physi- 
cal processes and the free parameters in the models limits 
the achievements. 

There is a class of galax y groups, dubbed as fossil 
groups (jPonman et al.l ll994!). which are the archetypal 
relaxed systems and arguably the end product of the 
galaxy mergers within the gr oup. The s e lection cri- 
teria for fossils is outlined by iJones et"al] (|2003D . i.e. 
groups with an X-ray luminosity of Lx^boi > 0.25 x 10*^ 
h~^ergs~^, and a minimum luminosity difference of 2 
magnitude between the first- and second-ranked galax- 
ies (AMi2 > 2), within half the projected i?200 radius, 
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where -R200 is the radius within which the mean density 
is 200 times the critical density of the Universe. 

Implementation of cosmological simulations in which 
the evolution of the halos can be traced, has given 
strong evidences that fossil groups, or galaxy sys- 
tems with large luminosity gap AA'/i2, on average 
form ear lier than galaxy groups with small luminos- 
ity gap HDariush et alj 120071 |2010[) . The early forma- 
tion epoch for fossil groups had been argued previ- 
ously based on observations of about a dozen such sys- 
tems. They included the study of X-ray scaling relations 
(jKhosroshahi et all 120071) as well as the morphological 
studies of the brightest group galaxies ( BGGs) in fossil 
groups (jKhosroshahi. Ponman fc Jones 1 120061 . 

A statistical comparison between the properties of ob- 
served fossil groups and those in the cosmological simu- 
lations requires dozens of such systems to be identified 
observationally which is currently unavailable due to the 
nature of existing X-ray surv eys (often sha l low o r very 
limited in angular coverage). iSantos et al.l (|2007[ ) have 
cross-correlated optical and the Sloan Digital Sky Survey 
(SDSS) and the ROSAT All Sky Survey (RASS) (Voges 
et al. 1999) and identified 34 fossil group candidates cov- 
ering a wide range of redshift. The applied methodology, 
however, does not produce a complete sample of fossil 
groups. 

Although conventionally an X-ray threshold has been 
applied to the IGM in fossil groups a lot can be learned 
from the selectio n based on opt ical criterion alone, e.g. 
"optical fossils" (jPariush et al.| [2007,). Identifying low- 
mass fossil systems is hampered by the fact that groups 
are under-represented in exis ting X-ray catalogs, how- 
ever, based on simulation data lDariush et all ()2007f ) have 
shown that even applying X-ray criteria, the fraction of 
late-formed systems that are spuriously identified as fos- 
sils is ~ 4 — 8 per cent, almost independent of halo mass. 

The main driver of this study is to find out the extend 
at which the luminosity gap in low mass groups can shed 
light on the formation of these systems and quantify the 
observed properties related to the luminosity gap which 
can be compared with semi-analytic models. As the large 
luminosity gap between the galaxies in a galaxy group 
can have a statistical origin, it would be therefore use- 
ful to find out to what extend the large luminosity gap 
is a representative of the evolutionary processes as op- 
posed to have been originated statistically. Furthermore 
we provide observational measures which can be used to 
constrain the semi-analytic models of galaxy formation. 

The paper is organized as follows; Section 2 describes 
the databases and the simulations. Algorithm used to 
extract the sample of groups in observation is described 
in section 3. Biases and completeness of the sample is 
discussed in section 4. In section 5 we present the derived 
group properties. Selection criteria for identifying fossil 
group is discussed in section 6. Results are presented 
and discussed in sections 7 to 9 with concluding remarks 
presented in section 10. 

For this study we assume, h = 0.7, flm — 0.3, and 
= 0.7. 



2. DATA 
2.1. Observation 



We use the legac y archive of the lates t data release 
of the SDSS, DR7 (iAbazaiian et al.ll2009D which covers 
8,423 square degrees in imaging data (which contains 
roughly 360 million distinct photometric objects) and 
the spectroscopic survey mapped 8032 square degrees 
(1,640,000 objects with measured spectra). We confine 
the study to those objects identified in SDSS as galax- 
ies. In this study, we use both the photometric and the 
spectroscopic data. 

Initially we restrict the study to galaxies within the 
redshift range of 0.02 < z < 0.17, retaining only those 
with clean spectra (fiag=0) providing us with a sample 
of 601,981 galaxies. The lower limit on the redshift is 
dictated recalling that, objects with z < 0.02 are often 
stars and the upper limit on redshift is placed in order to 
confine the study to the range where the spectroscopic 
sample is complete. It worths stressing that, the main 
galaxy spectroscopic sample of SDSS is restricted to Pet- 
rosian magnitude, rpetro < 17.77. 

We extract objects with rpetro < 21.0, identified as 
galaxy, from the photometric sample of SDSS DR7. For 
galaxies fainter than rpetro = 21. 0, the err or in ap- 
parent magnitude increases rapidly (jOvaizu et a l. 2008) . 
For objects satisfying the above conditions, we retrieved 
the P etrosian magnitudes (jPetrosian II1976I : iStrauss et al.l 
I2OO2D and K-corrected th em using the method described 
in lBlanton et al.l ()2003aD . 

2.2. Simulation 

In this paper we use th e galaxy group cata log from 
the Millennium simulation (jSpringel et al.lf2005D and the 
associated s e mi-an alytic model of galaxy formation by 
IBower et al.l (|2006f) . Below we briefly describe these sim- 
ulations. 

The Millennium simulation utilizes a ACDM cosmo- 
logical model to follow structure formation from z — 127 
up to the present epoch. Starting with an inflation- 
ary, dark-matter dominated universe, structures form 
through the bottom- up hierarchy, i.e., starting with small 
scale density fluctuations resulting in large scale struc- 
tures we observe today. Evolution of 2160'^ particles each 
with 8.6 X lO^h-'^MQ has been followed fr om z — 127 
up to present day (64 time-slices of positions and ve- 
locities were stored, separated logarithmically between 
z = 127 and z = 0) through a comoving box of 500Mpc 
on each side, (iSpringel et al.l I2005D with fix = 0.75, 
flM = 0.25, = 0.045, h = 0. 73, n = 1 and as = 0.9 
based on WMAP observations, ('Spergel et al."2003') and 
2dF galaxy redshift survey, (|Colless et al.„ .2001') . The 
evolution of a structure has been followed only if it 
composed of at least, 20 particles (equivalent mass is 
1.72 X 10^°h-^MQ) at that epoch (jSpringel et al.| [2005V 
In addition, friends-of-friends algorithm was exploited to 
extract haloes with densities at least 200 times the criti- 
cal density and substructures were identified u sing S UB- 
FIND algorithm developed bv ISpringel et al.l (|2001f ). 

The underlying dark matter haloes of the Millennium 
simulation was used to simulate the growth of galaxies, 
by self-consistently implementing a semi-analytic model 
of ga laxies on the outp uts of the Millennium simula- 
tion (jBower et al.ll200"6l ). This semi-analytic model ac- 
counts for the feedback from supernova explosions as 
well as AGNs in modeling the massive halos. In addition, 
AGN feedback has been considered as an operation which 
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quenches the star-formation and star-bursts are triggered 
both by merging and disk-instabilities. iBower et alj 
(|2006D predicts the luminosity function of galaxies in B 
and K band for the present epoch as well as the present 
day color distribution. For the purpose of this study, 
we use the catalogue of lBower et al.l ()2006l ) semi-analytic 
model and retrieved the absolute r-band magnitude for 
group members, the halo mass of the group as well as the 
mass of the subhalo containing the BGG at the present 
epoch. 

2.3. Monte Carlo Simulation 

The idea behind performing the Monte-Carlo simula- 
tion is that the luminosity gap between the two brightest 
members (or between any other members) of a galaxy 
system can also have a statistical origin. 

We expect large luminosity gaps to appear preferen- 
tially in groups with few members. In order to under- 
stand the extend at which the luminosity gap with statis- 
tical origin affects our analysis and results, we perform 
a Monte-Carlo simulation by drawing random sample s 
from an underlying Schechter function (jSchechter II1976I ) . 

We generate a sample of ^ 10^ galaxy groups with 
different number of members within the completeness of 
the observed sample described in section 4. We choose 
the underlying luminosity function to have Schechter 
parameters of a = — 1.08 and M* = —21.62 as in 
iZandivarez et al.l ()2006( ) , with galaxy luminosity ranging 
from Mr = —16 to —26. The chosen values are the best 
fit Luminosity function for groups based on the r-band 
SDSS data. 

3. SELECTING GROUPS IN OBSERVATIONAL 
DATA 

3.1. The Algorithm 

To identify galaxy groups, we use the m odified friends- 
of-frie nds (mFoF) algorithm develope d by iParedes et ah 
j[l99l) which is a modification of FoF ()Huchra fc Gellerl 
[1982:), one of the most frequently applied methods for 
finding structures in redshift surveys. The starting point 
in this algorithm is that every galaxy can be the center of 
a group of galaxies. The algorithm, however, will gradu- 
ally lead to the most probable center, for a given search 
radius. 

Given this, we start by looking for companions of each 
galaxy within a specified search physical radius, r^. A 
galaxy is companion to the chosen galaxy if it falls within 
the search radius in the plane of the sky and within the 
redshift interval of Az. For this study we choose rg to 
be the equivalent of 500 kpc at the redshift of the galaxy 
in the spectroscopic sample and Az=0.002. 

For the first generation of the groups, we count the 
number of companions with the above constraints for 
each galaxy in the spectroscopic sample described in sec- 
tion 2.1. Naturally most of the members will be shared 
by galaxy groups found in the first generation. Com- 
paring overlapping groups, those with fewer number of 
member will be removed from the list if they fall within 
the measured from the center of the richer group and 
Az=0.002. This results in the second generation of the 
galaxy groups. Having this, we calculate the mean red- 
shift and luminosity weighted center of each of the surviv- 
ing groups and the process of removing the small groups 



overlapping with larger groups continues until there is no 
common members between any two groups of galaxies. 

As it was described, mFoF works with two free 
parameters, rg and Az, and their valu es are set depending 
on t he scien ce goals of the s tudy, e.g. iMiller et al.l ()2005[ ) 
and Ivon der Linden et al.l ()2007l ) have been used r^ = 
IMpc ; Az=0.002 and = 2Mpc ; Az =0.01 to iden- 
tify cl uster of galaxies, respectively. Also, lEinasto et al.l 
(|2005D adopted = 0.5Mpc ; Az=0.002 to make cata- 
logue of group/cluster. Here we aim at studying galaxy 
groups and low mass haloes for which we find our choice 
suitable in comparison to the existing galaxy group cat- 
alogs. We extracted 81614 and 6956 groups with more 
than 2 and 5 spectroscopic members (richness), respec- 
tively. 

There exist several catalogues of groups and clusters 
of galaxies extracted from the SDSS data: Catalogues 
by Goto (2005, DR2), Miller et al. (2005, DR2), Mer- 
chan and Zandivarez (2005, DR3), Berlind et al.(2006, 
DR3), Zandivarez et al. (2006, DR4), Yoon et al. (2008, 
DR5) and Tago et al. (2008, DR5). In each of the 
above mentioned catalogs, criteria implemented through 
the group-finder algorithm has been chosen on the ba- 
sis of the type of the corresponding study. To check if 
our extracted groups, based on conditions adopted here, 
overlap with those from previous studies, we performed 
a cross-matching between our extracted catalog and two 
other popular group samples in the literature: the Abell 
cluster (JAbeUj 1958) and C4 (Mille r et all l200 5f) cata- 
logs. It worth stressing that although C4 is known as a 
"galaxy cluster catalog" , it has a wide mass range making 
it suitable for the comparison with our "galaxy group" 
catalog. 

For the red shift of Abell clusters we used published 
data by Strub ie~et al.l ([1987) in which there are 838 with 
determined redshifts out of 2712 Abell clusters. There 
are 318 Abell clusters that lie in the area of sky cover- 
age for the SDSS DR7 and have redshifts in the range of 
0.02 < z < 0.17. For our groups of richness > 5 (spec- 
troscopic members), 195 out of 318 Abell clusters (61%) 
are matched within 10" of our group center. In addition, 
we selected the groups of richness > 2 and in this case, 
283 out of 318 AbeU clusters are matched (88%). The 
C4 catalog was generated using a cluster finding algo- 
rithm which identifies clusters as over densiti es in a seven 
dimensional space of position and colors ([Miller et al.l 
2005). The catalog contains 748 clusters with richness 
more than 10 members brighter than r = 17.7. Applying 
the above search criteria we found 549 of 748 C4 clus- 
ters, 73% within 0.02 < z < 0.17. Reducing the number 
of galaxies per group, > 2, we would find 96% of C4 
clusters, which shows a very large overlap between our 
extracted groups and the C4 catalog. 

3.2. Color-Magnitude Relation 

As noted earlier, the main spectroscopic sample of the 
SDSS is complete to the Petrosian r < 17.77. How- 
ever, the photometric sample reaches lower luminosities. 
Moreover, the tiling algorithm use d in SDSS spectro- 
scopic survey ([Blanton et al.ll2003b( ). leaves some galax- 
ies uno bserved spe c trosco pically because of fiber colli- 
sions. [Yoon et all ([2008D estimated the spectroscopic 
completeness of the SDSS DR5 to be fspec ~ 65% 
for rich clusters. In order to have a realistic estimate 
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Fig. 1. — Color-magnitude relation for members residing within a 
centered at the BGG position with a radius of 0.5 Mpc at the BGG 
redshift. The 2(t boarders and restrictions on apparent magnitude 
are also illustrated. 

of the optical luminosity and the richness of the ex- 
tracted groups we use the color magnitude relation (here- 
after CMR) to identify additional group member can- 
didates. Observational evidences have shown that the 
bulk of the early-type galaxies in clusters and groups 
lie along a linear CMR. This relation is shown to have a 
smaU scatter (e.g-.IBower et al.l ll992t [Peebles et alll2002l: 
IKodama et al.lllQQSl) and can be used as a robust method 
for finding galaxy systems (e.g., [Gladdcrs et al. 200&). 
Although there are uncertainties in the CMR method of 
group membership identification, we use the photometric 
members only to estimate the global properties, such as 
the number of galaxies per groups and the total luminos- 
ity of the group. 

In order to incorporate the photometric data into the 
analysis, we consider all the photometric and spectro- 
scopic galaxies within an angular distance equivalent to 
500 kpc around the group center if they fall on the CMR. 
It should be noted that, group center is defined as the 
brightest member of the group. For this purpose, we fit a 
line to the color-magnitude diagram of the spectroscopic 
members if the group has at least three members with 
Mr < —19.0. We determine the slope and intercept of 
the CMR for each group. We consider only those photo- 
metric (m^ < 21) and spectroscopic galaxies which fall 
within ±2(7 of the color magnitude relation defined by 
spectroscopic members. An example of a fitted CMR to 
color magnitude diagram for one of our groups is shown 
in Fig. m 

Member selection using the color-magnitude relation 
method is hampered by miss-identification of forground 
and background galaxies as group members. In princi- 
ple, miss-identification of members will affect the total 
luminosity and in turn the estimated mass (§5.2). In or- 
der to quantify the uncertainty in halo mass estimate we 
randomly selected field regions in the SDSS. We chose 
galaxies within a radius comparable to that of groups in 
each field region and applied color-magnitude selection. 
The contribution of the field galaxies on the halo mass 
of groups was found to be negligible. 

We assume all the photometric group members found 
using the CMR to have the mean redshift of the group 
defined by the spectroscopic members. 
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Fig. 2. — 4113 groups with M^^bgc !5 ^22 are shown in red- 
shift vs. Ami2 diagram. The green rectangle area represents the 
selected region which is free of biases. 



4. BIAS CORRECTION AND COMPLETENESS 

As noted in §2.1, we have used the main spectroscopic 
sample and photometric legacy survey of SDSS DR7 
as our spectroscopic and photometric input. Statistical 
properties of the input data is discussed here along with 
the constrains resulted from the completeness and bias 
removal from th e samp les. For the photometric sample, 
lAbazaiian et al.l (|2Q09D demonstrated that the sample is 
95% (95% completeness limit for point sources) complete 
to r « 22.2. Whereas, the main spectroscopic sample 
of SD SS DR7 is comp lete to a r-band petrosian magni- 
tude (IPetrosian 1119761) of r :^ 17.77 , corrected for Galac- 
tic extinction (ISchlegel et al.l 119981: IStrauss et "all 120021: 
lAbazaiian et al.ll2009D . 

The limiting magnitude of r ~ 17.77 is equal to — 
-17 aX z = 0.02 and Af,. ~ -21.8 at z = 0.17. As we 
are interested in large luminosity gap, i.e. Ami2 > 2 
magnitude, the above completeness limit will naturally 
bias our analysis toward low luminosity gap at higher 
end of the redshift distribution. This is shown in Figl^l 

To overcome this, we apply the Mr< —22 constraint 
to the BGGs both in the observed SDSS groups (4113 
groups) and in the simulations (18843 groups). The con- 
strain on the BGG luminosity is also consistent with the 
value adapted for the M* parameter in Schechter lumi- 
nosity function in r-band. In addition this selection helps 
to compensate the absence of X-ray data, given the corre- 
latio n of the BGG luminosity with the IGM X-ray emis- 
sion (|Khosroshahi et all [20071: lEllis fc O'SuiirTOl1 [2006V 
and thus reduces contamination from individual galaxies 
that are not associated with a group scale halo. 

To ensure that a luminosity gap up to 2.5 mag can 
be observed without a bias, the second brightest galaxy 
should be at least as luminous as -19.5 mag which corre- 
sponds to the limiting apparent magnitude of the spec- 
troscopic sample, i.e r ~ 17.77 mag at z = 0.065. 
In other words, our luminosity gap estimates less than 
2.5 magnitude should be bias free for groups within 
0.02 < z < 0.065. We therefore restrict the observed 
galaxy sample to the above redshift range. 

Applying the above mentioned restrictions (Ami2< 
2.5 and z < 0.065) resulted in 109 fossils out of 620 
groups (see. §6). 
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ESTIMATING PHYSICAL PARAMETERS FOR 
GROUPS 



galaxy group, we have derived a number of quantities 
both for simulated and observed groups which is dis- 
cussed here. 

5.1. Velocity Dispersion and Virial Radius 

According to iBeers et all (|1990f ). a reliable method 
for measuring the velocity dispersion in groups with 
few m embers is the so-called gapper estimator (e.g. 
lYang e t al. 2005). Since our groups catalogue contains 
low member groups we use this method to estimate the 
line-of-sight velocity dispersion of each individual group. 

The method involves ordering the set of recession ve- 
locities Vi of the N group members and defining gaps 
as 
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l,2,...,N- 1. 
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The rest-frame velocity dispersion of our groups is then 
given by 
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Fig. 3. — The distribution of the galaxy groups in the plane of 
total mass of the group and the BGG luminoisty. Total mass of 
the simulated groups are estimated based on total luminosity Eq 
6 (left) and dark matter particles (right). The mass of observed 
groups are based on the total luminosity. It is noticeable that, total 
mass obtained through total luminosity is in a better agreement 
with the observations. 
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where the weight is defined as Wi — i{N ~ i) ■ As the 
central galaxy in each group, is assumed to be at rest with 
respect to the dark matter halo, the estimated velocity 
dispersion has to b e corrected by a factor of y/N/{N — 1) 
(|Yang et al.ll2005l ). This resuhs in a final velocity disper- 
sion given by 



N -1 



As we will see in the next section, an estimation of 
the virial radius of the galaxy groups is also required. 
Using the Virial theorem, Girardi et al. (1998) derives an 
approximation for the virial radius of a spherical system 
as: 



Ry 



To 



Mpc 



where ay is the velocity dispersion of group members. 

On the other hand, for the simulated groups, given that 
the halo mass is derived directly by adding the mass of 
the dark matter particles, Rmr is calculated as follows: 



Ry 



GMhalo 



Mpc 



where M^aio is the total mass of the halo and ay de- 
notes the velocity dispersion of group members. 

5.2. Halo Mass 

For observed groups, the total luminosity was used 
as diagnostic of halo mass, as appeared to be a better 
mass estimator than the line-of-sight velocity dispersio n 
(IPonesso et al.l[2005l : iMilTer et al.ll2n05l: ILin et al.ll200l : 



log{h-^M) w -2AQ+lAUog{h''^Lr) 



where is the total r-band luminosity of the galaxy 
Wigi. (2) group, defined as the sum of the r- band luminosity of 
member galaxies with Mr < —19.5, ^Milosavlievic et aLl 
[2006h . 

In addition to the measurement of Mhaio in the sim- 
ulated groups, we also estimate the halo mass using the 
empirical relation (Eq. 6) for consistency. In Fig. |3l 
the two methods of mass estimation are compared. Left 
panel compares the distribution of the halos in the plane 
of total gravitational mass based on Eq. 6 and the lu- 
minosity of the BGG for both the simulated and the ob- 
(|3) served systems. The same is shown in the right panel 
with the mass of the simulated halos obtained from the 
dark matter only. As expected, a better consistency is 
achieved when the luminosity based estimation is applied 
to the observations and the simulations. Therefore this 
mass indicator is used in the our analysis to avoid any 
systematics. The same method was applied in assigning 
a mass to the halos in the Monte Carlo simulations (see. 
§2.3). 

5.3. Richness 

To obtain the number of galaxies per group we applied 
Mr < —19.5 criterion to the observed and the simulated 
groups (Millennium and Monte-carlo). The distribution 
of richness for the observed and the simulated groups is 
shown in Fig. [4l This plot illustrates that, richness of 
(|5^ simulated and observed groups are different which might 
be related to different group selection procedures in simu- 
lation and observation data. The group selection in sim- 
ulations is based on a full knowledge of the three dimen- 
sion distribution of dark matter subhalos (e.g. galaxies), 
while as noted earlier, the group selection in observa- 
tions is based on two dimension distribution of galaxies 
complemented by slices the redshift space. However, the 
difference in shape of the distributions has no impact on 
our conclusions. 

In order to study the impact of biases on statistics of 
number of galaxies per group, generated groups through 
(6) the Monte Carlo simulation had to span a large range 
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Fig. 4. — The normalized histogram of number of galaxies per 
group is compared between simulated and observed groups. The 
K-S test shows that the probability of two samples selected from a 
common distribution is ~ 10^^^. 

of richnesses (10* simulations were carried out for each 
richness class of group). 

6. SELECTION OF FOSSIL GROUPS 

In low-density environments, the merging of compact 
groups can lead to the formation of the fossil groups 
(jVikhlinin eFaH 119991: Fjones et alJ [20031. These class 
of systems have masses which are comparable to those 
of normal groups and clusters of galaxies. Observation- 
ally, the identiiication of fossils have so far been based 
on the definition of fossil group suggested by Jon es et al.l 
(|2003f) . In iJones et all ()2003l ) definition, a fossil group 
has an X-ray luminosity of Lx,boi > 0.25 x 10*^ 
erg s~^, and the dominant galaxy is at least 2 magni- 
tudes brighter (in r-band) than the second ranked galaxy 
within 0.5i?2oo of the group. Assuming a circular orbit, 
the L* galaxies merge with the central galaxy within a 
few Gyr due to orbital dec ay as a result of the dynamical 
friction () Jones et al.ll2003D . 

The existing X-ray surveys overlapping the area cov- 
ered by the SDSS surv ey are either wide and shallow 
(ROSAT All Sky Survev. lt^s et alJlTool or very small 
and deep. As a result, facing the lack of suitable X-ray 
data we limit the study to opt ical fossils or galaxy groups 
with a large luminosity gap (iDariu sh et al.l [20071 ) . The 
X-ray criterion is conventionally used to ensure that only 
group-mass halos have been selected. As noted earlier we 
compensate this by the introduction of lower limit in the 
BGG luminosity. 

The luminosity gap was measured within half a virial 
radius, Rvir, derived using equations 4 a nd 5. This 
is slightly different to the definition of the IJones et alJ 
(|2003). however the difference applies to both the simu- 
lations and observations, consistently. 

Similar to observational constrains discussed in §4 ,the 
same conditions was applied to simulation data to select 
optical fossil groups with the excepti on of the X-ray l umi- 
nosity criterion. Also following Dari ush et al.l ()2007D . we 
limit the analysis to groups with log{Mfiaio/ Mq) > 
13. 

These conditions resulted in 1993 fossil systems out 
of 16688 groups, or roughly 12%. Furthermore, fossil 
groups through the Monte Carlo simulation, were defined 
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Fig. 5. — The distribution of r-band absolute magnitude of fossil 
BGGs is compared to the same in control groups for observational 
and simulation samples. Fossil BGGs are brighter than control 
BGGs. The K-S test shows that the probability of control and 
fossil samples being drawn from similar distribution are 1.2 X 10~^ 
and 3.6x 10~®* in observation and simulation samples, respectively. 
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Fig. 6. — The distribution of BGG absolute magnitude over group 
halo mass is compared between fossils and controls for both simu- 
lated and observed groups. 



as groups Atoi2 > 2 and Mr^BGG < —22 in accordance 
with definitions for observation and simulation data. 

In order to compare the properties of fossil with non- 
fossil groups, we introduced control groups with Ami2 < 
0.5 and satisfying the remaining conditions mentioned 
above (e.g. Mr.BGG < —22), in observations and simula- 
tions. This provides us with 85 and 2419 control groups 
in observed and simulated catalogs, respectively. 

7. FOSSIL GROUPS: BGG LUMINOSITY AND 
HALO PROPERTIES 

As noted earlier, fossil groups are distinguished by a 
large gap in the luminosity of the two brightest members 
indicating the absence of L* galaxies. If the BGGs in fos- 
sil groups are the product of the mergers of L* galaxies, 
then they are expected to be statistically brighter than 
the BGGs in non-fossil groups. 

Fig. [5l shows the distribution of the r-band absolute 
magnitude, M^, for fossil and control BGGs. As ex- 
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pected, the fossil BGGs are brighter than control BGGs 
in both the simulated and observed group s. This is also 
shown in a study by iDiaz-Gimenez et alj (|2008[ ) where 
they use simulations to demonstrate that the fossil BGGs 
on average accrete larger stellar mass as a result of earlier 
formation epoch and larger frequency of mergers. 

In Fig. ini we compare fossil and control groups in the 
plane of group mass and the BGG luminosity. It illus- 
trates that, for a given r-band luminosity of the BGG, 
control group halos are statistically more massive than 
fossil group halos. 

8. LUMINOSITY GAP STATISTICS 

As noted in §2.3, for some of the properties studied 
in this paper, there exists a probability that, random 
processes contribute by part to the observed trends. In 
the following, we will specifically study the contribution 
of random processes on halo properties of groups. 

8.1. The abundance of fossil groups 

In Fig. [7]-(a), we show the probability distributions of 
fossil groups for a given BGG r-band luminosity, in obser- 
vations as well as in the Millennium simulation and ran- 
domly generated groups through the Monte Carlo simu- 
lation. 

As shown in Figl7]-(a), within the uncertainty, the 
probability of finding large luminosity gap (fossil groups) 
is nearly independent of the luminosity of the brightest 
galaxy, in both the observed and simulated galaxy sys- 
tems. However the same increases with the luminosity of 
the brightest member in Monte-Carlo generated groups. 
The shape of the luminosity function is what drives the 
trend seen in the Monte-Carlo generated systems. As the 
probability of finding low luminosity galaxies is larger in 
comparison to luminous galaxies in the Schechter func- 
tion, thus, for a fixed M* more luminous BGGs would re- 
sult in larger luminosity gap, statistically. The fact that 
such a trend is not pronounced in the observed sample 
and in the cosmological simulation, supports the argu- 
ment that the large luminosity gap is unlikely to have a 
statistical origin. Moreover, as shown in the Figure [T]- 
(a), the fraction of fossil groups in the observations and 
the simulation is larger than the same predicted by the 
Monte-Carlo simulations. 

Fig. [7]-(b) shows the same as a function of group rich- 
ness. As expected, the probability of finding large lumi- 
nosity gap (Ami2> 2) decreases with increasing richness. 
However the trend is the fastest for the Monte-Carlo gen- 
erated groups while observed groups show a weaker de- 
pendency on richness. 

In the regime of low richness, the probability distribu- 
tions of the observed, simulated and randomly generated 
groups coincide. As a result, it appears that the luminos- 
ity gap in such systems is driven by random processes, 
i.e., evolutionary mechanisms have played negligible role 
in forming the luminosity gap. This is supported by the 
hierarchical structure formation paradigm in which the 
low mass hal os (poor galaxy groups ) are recently formed, 
statistically. iDariush et al.l (|2007[ ) concluded the same 
based on a comparison of a simulated (Millennium) and 
Monte-Carlo generated groups. This study extends their 
findings to the observed groups. 

On the opposite, in richer groups, the Monte-Carlo 
prediction of the large luminosity gap incidence is sig- 



nificantly lower than those in cosmological simulations 
and the observations suggesting that the luminosity gap 
in richer systems is predominantly driven by evolution- 
ary processes. Moreover, it is noticeable that the cos- 
mological simulation based on semi-analytic model of 
I Bower et~al. (2006), predicts a lower probability for fos- 
sil groups, in the entire richness range, than the same 
in the observed sample. Furthermore, choice of group 
finding algorithms in observations and simulations could 
also contribution to this difference. The probability dis- 
tribution therefore helps to probe the accuracy of the 
semi-analytic models. 

In Fig. [7]-(c) we show the fraction of fossil groups in 
each mass bin. For Monte Carlo generated groups, it is 
prominent that, the fraction of fossil groups in each mass 
bin, decrease as the halo mass increases. As noted earlier, 
masses of each Monte Carlo generated groups were as- 
signed according to the total luminosity method. There- 
fore, we expect that, both the richness of the groups as 
well as the luminosity of bright members determine the 
halo mass. For fossil groups (and generally for systems 
with large luminosity gap), the BGG luminosity is a more 
defining factor because of the presence of the large lumi- 
nosity gap. Thus, fossil fraction in each mass bin, is 
strongly linked to the BGG luminosity and the richness 
of the group. 

Comparing the probability distributions of observed, 
simulated and randomly generated groups, illustrates 
that the three distributions coincide in low mass end. 
Consequently, low mass fossil groups are predominantly 
statistical i.e., evolutionary processes are not playing a 
significant role in forming the luminosity gap in these 
systems. On the contrary, at the high mass end, fraction 
of fossils are substantially higher than what is predicted 
by the Monte Carlo simulation both in the cosmolog- 
ical simulations and in the observations, indicating the 
role of evolutionary mechanisms in forming Ami2. Inter- 
estingly semi-analytic model adopted for this study pre- 
dicts a lower fraction for fossil groups across most of the 
mass/richness range. It is hard to point out at a single 
physical process employed in the simulations to interpret 
the offset between the fossil fraction in the observations 
and simulations, as this could be a consequence of vari- 
ous assumptions in handling non gravitational processes 
such as the starformation efficiency, AGN feedback and 
cooling. These could be coupled with gravitational pro- 
cesses such as the merger rate. 

We propose this to be an observational constraint for 
the semi-analytic models in addition to the luminosity 
function and the color of galaxies. 

8.2. Luminosity Gap Probability Distribution 

We study the probability distribution of luminosity 
gap. The distributions are shown for two richness 
regimes, richness < 5 and richness > 15 in Fig. |S1 
It is noticeable that, although the probability distribu- 
tions in observed and simulated groups follow a similar 
trend (see below) , they deviate from Monte Carlo predic- 
tions at Arni2> 1.2 and Ami2> 0.75 for low and high 
richness respectively, such that the fraction of observed 
and simulated groups with a given Arni2 are larger than 
randomly generated groups. 

We interpret this as an impact of evolution: assuming 
that, the primordial distribution of Ami2 was what is 
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Fig. 7. — The upper panel represents the percentage of fossil 
groups in each absolute magnitude bin, for randomly generated 
as well as observed and simulated groups. The percentage of fos- 
sil groups in each richness bin is illustrated in the middle panel. 
The lower panel shows the fossil percentage in each halo mass bin, 
for simulated, observed and randomly generated groups. In richer 
and more massive groups the evolutionary process dominate the 
random processes. 



predicted by Monte Carlo simulation, the evolutionary 
mechanisms, always work in the direction of increasing 

t he luminosity gap. 

Ivon Benda-Beckmann et all (j2.008) suggested an in- 
verse process; transforming systems with large luminos- 
ity gaps to those with low Ami2- In this picture, there 
is always a probability that, new galaxies fall into groups 
and hence disrupt the luminosity gap. For instance, they 



Number of Galaxies per Group > 1 5 




Fig. 8. — Normalized histogram of Ami2 for two separate rich- 
ness regimes (see. §8.2). 

have shown that, luminosity gap is strongly affected by 
galaxy infall and hence is a transient characteristic of 
galactic systems. Other processes may also affect the lu- 
minosity gap, including mergers between galaxy groups. 
However when two galaxy systems merge, they end up 
in considerably more massive systems and hence with 
higher richness than the progenitors. This also works in 
the direction of decreasing the luminosity gap. 

Fig. IS] shows the relative efficiency of evolutionary pro- 
cesses resulting in larger luminosity gap compared to in- 
verse processes which works in the opposite direction. 
The difference between the two panels in this figure also 
indicates that poor systems (low richness) on average are 
those that are closer to the Monte-Carlo predictions and 
therefore could have recently formed and least affected 
by evolutionary processes. Another possible explanation 
would be that, in such systems the evolutionary processes 
simply operate more slowly. 

9. NEW CRITERION FOR FINDING FOSSIL 
GROUP; Amu 

In a recent study based on the Millennium simula- 
tion bv iDariush et al.l (|201Cll ). a new indicator based on 
gap magnitude between first and forth brightest galaxies 
within half of the virial radius of the group Amu ^ 2.5 
is presented to be a more efficient probe of identifying 
early-formed halos than conventional definition of fossil 
g alaxy group, i.e . Ato i2 > 2. 

IDariush et al.l (|201C1( ) find that the mass assembly his- 
tories of the halos identify by the two methods, on aver- 
age, are similar. About 90% of fossil groups which were 
identified according to Ami2 and Amu criteria in earlier 
epochs become non-fossils after ~ AGyr and the fossils 
phase itself lasts ~ IGyr. The main difference between 
the two methods se em to be in the effi ciency of finding 
early formed halos (|Dariush et al.ll2010l ). 

In order to evaluate the abundance of fossil groups ac- 
cording to these definitions, we compared the luminos- 
ity gap distribution based on simulation and observation 
data (our catalog). 

To be able to estimate the luminosity gap parameter 
Amu, groups with the following properties are selected. 
The extracted groups have at least four members within 
the half of the virial radius of the galaxy group and the 
brightest and fourth brightest galaxy should be at least 
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as luminous as —22 and —19, respectively. Applying 
above criteria, 163 groups were identified in our catalog, 
out of which 41 groups meet the Atoi4 > 2.5 criterion. 

Among these observed fossil groups (Ami4 > 2.5), 12% 
satisfy the conventional definition of fossil (Atoi2 > 2). 
Conversely, 56% of Ami2 > 2 satisfy the Ami4 > 2.5 
criterion. Hence, a large proportion of the population 
of groups identified with Ami4 > 2.5 criterion are not 
in common with those identified by Ami 2 > 2 . This 
is consistent with the findings of IDariush et al.l (|2010( ) 
based on simulations. 

We applied the same criteria to extract galaxy group 
in simulation data. Out of 16688 groups, 1993 and 2973 
groups satisfy the Ami2 > 2 and Ami4 > 2.5 groups 
respectively. 

In Figini we plot the R-band luminosity gap distribu- 
tion Atoi2 and Ami4 from our observations and simula- 
tions. Results show that the luminosity gap from simu- 
lation is in a fair agreement with observations and both 
present a similar distributions for the Ami2 and Atoi4 
luminosity gap. The fraction of groups with Ami2 > 2 
in observation and simulation are 17% and 12%, respec- 
tively while these values for Ami4 > 2.5 indicator are 
25% and 18%, respectively. 

The probability of finding fossil system according to 
these definitions in c urrent analysis is diff erent to the 
statistics reported by 'Dariush et al.|_( |2010[ ) which uses 
the group catalog of Yang et al. (20 071) ba sed on SDSS 
DRA. For instance iDariush et all J2010|) reports 2.0% 
and 2.1% in observation and simulation, respectively 
when using Atoi2 > 2. For Atoi4 > 2.5, they find 6.2% 
and 5.1%, respectively. The difference is due to the com- 
pleteness considerations in sample selection adapted by 
us. 

The most imp ortant difference be tween our group cat- 
alog and that of lYang et al.l ()2007l ) is the addition con- 
strain we have imposed on BGG luminosity and the lu- 
minosity gap (see. §4). Furthermore their groups are 
defined as systems which at least four members while 
we require the group to have two spectroscopic member 
within half of the virial radius. 

In §8.1, it was found that the probability of finding 
early-formed system, i.e. fossil groups Ami2 > 2, is 
higher in low mass systems. In FigfTUl the abundance of 
fossil groups base on conventional definition Ar7ii2 > 2 
and the new indicator Atoi4 > 2.5 are given as func- 
tion of halo mass. This shows that, in comparison to 
conventional fossils (Ami2 > 2), the fossil groups based 
on Arni4 > 2.5 are more abundant in low mass range 
log{M /h^^Mo) < 14.3. Above this mass limit, the 
statistics is poor and therefore hard to interpret. 

10. DISCUSSION AND CONCLUSION 

In this study we extend the earlier observational stud- 
ies of luminosity gap to low mass regime and provide ob- 
servational constrains for semi-analytic models based on 
this observable. We extract 620 groups in SDSS DR7 
utilizing mFoF algorithm and measure the luminosity 
gap between the two brightest galaxies in the group, 
Ami2. This results in 109 groups with Ami2 > 2, within 
half a vir ial radius, known a. s optical fossils groups. In 
addition, iBower et al.l (|2006f) semi-analytic model em- 
ployed on Millennium simulation was used to select 16688 
groups with 1993 fossil groups in the same manner as 
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Fig. 9. — Distribution of Ami2{a) and Am,i4(b) for observa- 
tional and simulation samples. 
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Fig. 10. — Fraction of fossil groups according to both definitions 
Ami2 > 2 and Ami4 > 2.5, as a function of halo mass, in obser- 
vational sample. 

applied to the observational data. A Monte carlo simu- 
lation was also performed to generate a large sample of 
luminosity gap statistics purely drawn at random from 
Schechter luminosity function to investigate the impor- 
tance of the random processes. 

We show that fossil BGGs both in observations and 
simulations are on average brighter than those in control 
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galaxy groups. 

Given a relatively poor constrain on galaxy velocity 
dispersion in groups (e.g. dynamical mass estimation) 
we adapted the total luminosity method for the estima- 
tion of the total gravitational mass of the groups based 
on empirical relations. We show fossil groups halos are 
found to be statistically less massive than those of the 
control sample for a given BGG luminosity. 

We confirm that the luminosity gap in systems with 
low richness is predominantly driven by random pro- 
cesses while evolutionary mechanisms are responsible for 
large luminosity gap in rich groups. 

We find indications that the fraction of fossil groups 
based on Ami4 > 2.5 criterion is higher compared to the 
same when adopting the conventional Ami2 > 2 crite- 
rion. This is valid mostly to t he low mass end. 

Following a recent study bv ISmith et al.l ()2010( ). base 
on a sample of massive galaxy clusters (~ IO^^Mq), we 
propose a test based on the luminosity gap in galaxy 
group, the butterfly diagram vs. Ami2, which is 
found to be able to discriminate between different galaxy 
formation models in groups and clusters when a large 
sample of groups with known luminosity gap is available. 
Extending their study to low mass groups, we compare 
620 observed groups and 16688 simulated groups in the 
plane of Mr vs. Ami2, Fig. [TT] The absolute magnitude 
of the brightest and second brightest galaxy within the 
group are plotted as a function of Ami2 for different halo 
masses. A linear regression is also presented for both the 
brightest and the second brightest galaxies as a function 
of Ami2. In addition. Fig. [TT] shows the variation of the 
slope and the intercept of the linear fit as a function of 
the halo mass for BGGs. 

We find that th e slope of the BGG brightness in 
I Bower et al.l (|2006f ) simulation is steeper than the ob- 
served BGGs, particulary in high mass regime. This 
i mplies that the AG N feedback in BGG is too weak 
in [ Bower et al.l (|2006l ) semi-analytic model (|Smith et al.l 
|2010( ). Further improvements have been made in the 
semi-analytic model recently (Bower et al 2008) driven 
by constraints from the IGM properties of the galaxy 
systems, which also required a re-tuning of the galaxy 
formation parameters in their earlier model to be able to 
reconstruct the galaxy luminosity function in the local 
universe. 

The discussion on the success of their changes in the 
context of the luminosity gap requir es a direct cq r apari- 
son o f the two semi-analytic models iBower et al.l (|2006l 
I2008| ). We find that the BGG brightness of observed 
groups, increases with halo mass. In contrast this trend 
is not clear in simulated groups. Furthermore, in high 
mass bin, the absolute magnitude of simulated BGG span 
^ 2.5 mag, in cont rast to the observed range of '-^ 1 mag. 
ISmith et all (I20I0I) argue that the large spread in the ab- 
solute magnitude of simulated BGGs is driven by higher 
efficiency of the conversion of the cold gas into stars. We 
note that in low mass bin, where the random processes 
are dominated, the spread in the luminosity of the ob- 
served and simulated BGGs are nearly the same (^ 0.5 
mag). 

The correlation between the slope of the line fit- 
ted to BGG in Fig. [11] and halo mass points to 
an increasing contribution of the evolutionary mech- 
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Fig. 11. — Upper Figure shows the distribution of r-band absolute 
magnitude for the first (open circles) and second (filled circles) 
brightest group member as a function of Ami2 for different halo 
masses. Figure at the bottom shows, the evolution of slopes and 
intercepts of brightest group member as a function of halo mass. 



anisms in forming the luminosity gap in more mas- 
sive halos. Giving the consider able space density of 
large l uminosity gap groups (e.g. Vikhlinin et al. '1999(, 
[Jones et al. 2003^. van den Bosch 2007. D 'O nghia et al] 
2005, von Be nda-Beckmann et al. I l2008inDariiIsh' et al.l 
i2007i ). it seems that the butterfly diagram is a simple 
way of testing the accuracy of the semi-analytic models. 
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